function y = source(x)
    y = pi^2 * sin(pi * x);
end